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The separability of the continuous- variable EPR state can be tested with Hanbury-Brown 
and Twiss type interference. The second-order visibility of such interference can provide 
an experimental test of entanglement. It is shown that time-resolved interference leads to 
the Hong, Ou and Mandel deep, that provides a signature of quantum non-separability 
for pure and mixed EPR states. A Hanbury-Brown and Twiss type witness operator can 
be constructed to test the quantum nature of the EPR entanglement. 
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1. Introduction 

The famous debate between Einstein, Podolsky, Rosen (EPR) and Bohr 1,2 , about 
the nature of quantum correlations of a bi-partite state has played a key role in 
the investigation of entanglement properties of light in modern Quantum Optics, 
and has initiated a new branch of physics called Quantum Information. EPR used 
in their arguments a wave function that exhibits a perfect correlations between 
positions and momenta of two massive particles (labeled a and b). In the position 
and momentum representations the EPR state takes the following forms 

^{x a ,X b ) ~ 6(x a - X b ), ^(p a ,p b ) ~ S(p a +p b ). (1) 

Quantum correlations of the entangled EPR state have been implemented ex- 
perimentally for massless particles: photons. In the case of light, the quantum me- 
chanical position and momentum observables are played by the electric field: phase 
and amplitude quadratures. 

The well known two-mode Gaussian squeezed states of light are physical real- 
ization of the EPR state. Those states lie at the heart of quantum cryptography ? , 
quantum information 4 and quantum teleportation 5 . 
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There are few efficient ways of producing Gaussian EPR correlated states. One of 
them uses a Kerr nonlinear medium in optical fibers to entangle phase and amplitude 
quadratures 6 . However, the nonlinearity has to be relatively small to achieve a 
Gaussian state. The other method uses a beam splitter as a nonlocal operation that 
creates entanglement. A typical setup consists of two initially separable amplitude 
squeezed beams which interfere at the 50/50 beam splitter. 

Recent applications of entangled two-mode squeezed states of light in quantum 
information processing have generated a lot of interest in mixed Gaussian states. 
Although there are mathematical criterions for the entangled properties for mixed 
two-partite Gaussian systems 7 , they are not easy experimentally realizable. 

It is the purpose of this paper to show that Hanbury-Brown and Twiss (HBT) 
interference can reveal the quantum nature of entanglement of two-particle mixed 
Gaussian states. We show that it is possible to construct a quantum witness opera- 
tor, closely related to the HBT interference 8 , that probe entanglement of a mixed 
EPR Gaussian state. Using Hong, Ou and Mandel (HOM) interferometry we dis- 
cuss the entangled properties of time resolved interference of the EPR state. We 
show that the deep in HOM interference can provide a useful tool to study quantum 
separability of continuous- variable Gaussian EPR states. 



2. Mixed EPR states 

A non-degenerated optical parametric amplification, involving two modes of the 
radiation field, provides a physical realization of the EPR state (1). The quantum 
state generated in this process has the following form 

oo 

l*> = S VP^\n,n), (2) 

n=0 

where p n — n n /(l + n) n+1 is a thermal distribution with a mean number n of 
photons in each mode. 

Using field quadratures eigenstates: |a; ,a;j,), we obtain that the wave function 
of such a system (2) is Gaussian and has the form 

(x a , x b \y) = le-^HX^.+^+V"!^")^ 1 ' . (3) 

In the limit of n — > oo, the two-mode squeezed state (3) becomes the original EPR 
state (1). The state (3) is not entangled only if n = 0. 

The simplest mixed generalization of the EPR state involves a Gaussian density 
operator p a b, being a Gaussian operator of the field modes described by the annihi- 
lation and creation operators (a, a*) and (6,6'). This Gaussian density operator of 
the two modes is fully characterized by its second moment expectation values of the 
modes. In the case of a mixed EPR state the only non-vanishing field correlations 
are 



{a' a) = (b^b) = n, (ab) = —m c , 



(4) 
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where n is a mean number of photons in each mode, and m c correspond to the 
amount of correlation between the two modes. For |m c | = y/n{fi + 1), the EPR 
state reduces to the pure state given by (3). 

As it has been discussed in several papers (see the tutorial 7 and references 
therein), this mixed EPR state is separable for n > \m c \ and its density operator 
can be expressed in a sum of product states (Werner separability criterion) 

p 3 b = ^Pip i {a)® p l (b), (5) 

i 

where p l {a), p l (b) are the density operators of the two modes and X^ft = 1> wrtn 
< Pi < 1. The mixed EPR state is entangled and non-separable if 



n 



< \m c \ < y/n(n + l). (6) 



3. Hanbury-Brown and Twiss Interference with EPR Pairs 

Second and higher orders of coherence of a light beam can reveal its quantum 
nature, which cannot be observed in Young-type interference experiments. Second- 
order interference involving intensity-intensity correlations have been first applied 
by Hanbury-Brown and Twiss in stellar interferometry 9 . In modern Quantum Op- 
tics, second-order interference has been used as powerful tool to study nonclassical 
properties of light 10 . 

Hanbury-Brown and Twiss (HBT) interference measures a second-order nor- 
mally ordered intensity-intensity correlation function. In Fig. 1 we have depicted a 
setup for HBT interference that involves two light beams with annihilation oper- 
ators a and b interfering at the beam splitter. A correlation between clicks of two 
detectors corresponds to a normally ordered intensity-intensity correlation function. 
The HBT interference exhibits a typical pattern of the form: 1 + cos((ySi — 1P2), 
where ip± and ip2 are phase differences between the beams a and 6 in front of the 
detectors. These phases include geometrical phases and phases due to the possible 
action of the beam splitter. In this setup the two phases can be controlled ex- 
perimentally. i/ 2 ) is a second-order interference visibility. For a classical source of 
light this visibility is always bounded: < 5. This classical limit is violated for 
single photons, as it has been shown in the pioneering experiments performed by 
Mandel 11 . 

We shall apply the HBT setup to study the mixed Gaussian EPR state. At the 
detectors (i — 1,2) the positive-frequency part of electric field corresponding to 
modes c and d (normalized to the number of photons) is as follows 

E^)= a -^f. (7) 
The corresponding field intensity operator at the screen is equal to 

I{ipi) = £(-)(<ft)£ (+) (<^) = ^(a Jf a + tfb±tfae- iVi ±aW Vi )- (8) 
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Fig. 1. Hanbury-Brown and Twiss second-order interference. 

From this expression, we obtain that the normally ordered second-order intensity 
correlation is 



1 



(:%>i) I{cp 2 ):} = - (: {at a + tfbf :) - 2{a s a tfb) cos(^i - <p 2 ) 



4 

+ (e^ 1 



^ 2 )(fet( a t a + tfb) a) + (e^ 1 - e l ^){a){a)a + tfb) b) 



(6t 2 a 2 )-e i ^+^)(at 2 6 2 ) 



(9) 



For the EPR beams in a mixed state, with mean n photons and correlation m c 
described by the relations from Eq.(4), the fourth-order field correlations needed in 
the HBT calculations are equal to 



(eta* a a) = 2n 2 , {tftfbb) = 2n 2 , (aW&) 



n 2 + \m c \ 2 . 



(10) 



As a result the HBT intensity correlation function (9), for a continuous variable 
mixed EPR state, has the following form 



(:Jfo>i) I(<p 2 ):) = \{'in 2 + \m c \ 2 ) {l - cos(^ - <^ 2 )} 
with the second-order fringe visibility equal to 

«( 2 ) 



n 2 + \m c \ 2 



3n 2 + \m c \ 2 ' 



(11) 



(12) 



For a separable (classical) output state we have: | < < i. In the case of no 
correlation in the output state |m c | = 0, the visibility is = -|, as it should be 
for a thermal state. For entangled states this visibility is quantum i.e., it violates 
the classical inequality: < \. 

In Fig. 2 and Fig. 3 we have depicted the visibility (12) as a function of the 
correlation parameter m c for two different values of n = 1,0.1. 

Based on the above visibility analysis one can introduce a HBT witness operator 

1 2aWab 



W (HBT) i . 

2 (:(/„ + J 6 ) 2 :)- 



The mean value of this operator, 

Tr{W^ BT W = 



2(3n 2 + K| 2 ) ' 



(14) 



I 



I 
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Fig. 2. The visibility (12) evaluated for fi = 1. The vertical grid line defines a border between 
entangled and separable EPR, states. For \m c \ > 1 the state is entangled and the visibility is 
greater than | - horizontal grid line. 
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Fig. 3. The visibility (12) evaluated for fi = 0.1. The vertical grid line defines a border between 
entangled and separable EPR states. For |m c | > 0.1 the state is entangled and the visibility is 
greater than | - horizontal grid line. 

is positive for separable mixed EPR states and negative for entangled mixed EPR 
states. 

4. Hong-Ou-Mandel Time-Resolved Interference 

The analysis of second-order interference, can be extended to the case of photon 
pulses which duration is long compared to the time resolution of photodetector. 
Such single-photon emitters are already in use 12 > 13 . 

A source produces single photon pulses with a relative delay which interfere at 
the beam splitter. Simple analysis of beam splitter's action on impinging photons 
shows, that both photons will leave the same outport of the beam splitter. A relative 
delay is small enough to keep the pulses partially overlapping. Even though the 
overlapping is not complete, its time duration is longer than the time resolution of 
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the detector. Therefore, the temporal effects cannot be neglected any more when 
detecting pulses and the time delay r between detected pulses at two separate 
detectors has to taken into account. Experiments investigating interference for such 
pulses with different frequencies have already been performed 14 ' 15 . 

We restrict our description of the single-photon wave packets in the space- 
time domain to one-mode fields with electric field operators: Ei + \t) — a(t)a and 
E { b +) (t) = 0(t) b. The phases ipi and if2 are included in the mode functions phases. 
These mode phases contain an additive stochastic component related to the random 
arrival time of the correlated photons at the beam splitter or detectors. The time 
dependent electric fields of the two-mode output state is described by the beam 
splitter transformation (see Eq. (7)) 



£ < + , (t) = ^'w+^'m ^ w _ M+'w-4> (15) 



The joint coincidence of detecting photons from mode c and d at delayed times 
t and t + t can be obtained directly from the following second-order temporal 
coherence function 

G {2) (t, t + r) = (t) 25<"> (t + r) (t + r) 25<+) (t)) . (16) 

Using the mode functions, with random phases in the output state, we note that the 
only contributing terms to the second-order coherence function are the following 

G (2 \t,t + r) = I{|a(t)| 2 | a (i + r)| 2 (aVaa) + |/?(t)| 2 |/3(i + r)\ 2 (tftfbb) 
+ [\a(t)mt + r)\ 2 + \a(t + r)\ 2 \m\ 2 

- (a(t + T)a*(t)(3(t)/3*{t + T) + c.c.) (o^+ofc)}. (17) 

As an example we shall use in the calculations a special model of phase fluctu- 
ations. We will assume that the mode functions have random phases typical for a 
stationary stochastic phase diffusion model. In this example, the only non- vanishing 
autocorrelations of the mode functions are 

(a*(t)a{t + T))=I a exp(-^j, (P*(t)f3(t + t)> = h exp (-^j . (18) 

In these formulas we have assumed that the statistical correlations are station- 
ary and Gaussian with coherence times r a and Tb- The stationary intensities are: 
(a*(t)a(t))=I a and (/?*(*)/?(*)> = h- 

Applying the above temporal correlations of the mode functions, combined with 
the EPR correlations of a mixed state given by Eq. (10), we obtain that the second- 
order coherence function is stationary (dependents only on time r) and has the 
form 

(r) = \ jn 2 {II + ID (n 2 + |m c | 2 ) IJ b exp (-^ - ^) } . (19) 
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We shall simplify this formula further assuming that the mode intensities are equal: 
I a = h = I- Setting + \ = \ the formula above reduces to 

T a T b T c 

G (2) (r) - y (3n 2 + K| 2 ) jl- W ( 2 ) e ~^j (20) 

with the second-order visibility equal to (12). From this formula we obtain that 
a joint coincidence probability to detect a photon at time t and another photon at 
time t + t is 

p(T) = l-v^ cxp(-T 2 ), (21) 

where T is a dimensionless time with r expressed in units of r c . 

This coincidence probability exhibits a typical Hong-Ou-Mandel deep 16 . In Fig. 

4, we have depicted the joint coincidence probability for a state with n = 1. The 
upper curve (n = 1 and \m c \ — 1) is a border between separable and nonseparable 
EPR states. The lower curve (n = 1 and \m c \ — \/2) describes HOM interference 
of an entangled pure EPR state. The region between the two curves corresponds to 
nonseparable mixed EPR sates. For n = 1 the minimum value the HOM deep for a 
pure state is achieved at a zero coincident rate and is equal to p m in — 0.4. 

For the EPR state given by Eq. (2), with a small value of fi, one can get a 
significant increase of the the deep. For a weak two-mode pure squeezed state (fi <C 
1) we can approximate the formula (2) by 

|*>~ V^|0,0) + V^|1,1). (22) 

We see that such a weak pure squeezed state is the same as the one discussed in 
the single photon experiment originally performed by HOM 16 . 

Fig. 5 illustrates the HOM deep for a state with n = 0.1. The upper curve 
(n = 0.1 and |m c | = 0.1) is a border between separable and nonseparable EPR 
states. The lower curve (n = 0.1 and |m c | = V0.11) describes HOM interference 
of an entangled pure EPR state approximated by the state exhibited above. The 
region between the two curves corresponds to nonseparable mixed EPR sates. For 
n = 0.1 the minimum value the HOM deep for a pure state is achieved at a zero 
coincident rate and is equal to p m in = 0.14. 

5. Conclusions 

In this paper we have shown that the Hanbury-Brown and Twiss type interference 
can be used as an experimental test of quantum separability for continuous variable 
EPR states. One can associate with this second-order interference a HBT witness 
operator quantifying entanglement for pure and mixed Gaussian EPR states. We 
have shown that time-resolved coincidence rate of joint photon detection exhibits the 
Hong, Ou and Mandel deep, that provides a signature of quantum nonseparability 
for pure and mixed EPR states. 
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Fig. 4. The plots of joint coincidence probability p(T) as a function of T for pure (the lower 
curve) and mixed separable (the upper curve) EPR states are depicted. The minimum value of 
p(T) for the pure state (n = 1 and \in c \ 2 = 2) is below i. For a mixed state from the border 
between entangled and separable states, (n = 1 and \m c \ 2 = 1), p(T) reaches the value of ^ but 
docs not exceed it. 
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Fig. 5. Similarly, as in Fig. 4, for a state with n = 0.1, probability p(T) exceeds the value of 
i for a pure state (|m c | 2 = 0.11, the lower curve). In this case the minimum value is smaller 
than for n = 1 and reaches 0.14. For the border state, between entangled and separable states 
(|m c | 2 = 0.01, the upper curve), the minimum value of p(T) is again equal to i. 
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